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Abstract 

The multi-body neutrino potential energy is analytically estimated for a 
spherical neutron star with a vector potential model for neutrinos. We show 
that the self-energy and the neutrino number of the neutron star coincide with 
the semi-classical values in the large volume limit, and confirm that there is 
no catastrophe in neutron stars with massless neutrinos. 
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1. Introduction 



The massless neutrino exchange leads to a long-range force and in particular to 
a multi-body potential force whose astrophysical effects were discussed long ago. p|, 
^ Recently Fischbach |^ argued that the multi-body potential gives an unphysically 
large self-energy to stellar objects like neutron stars and concluded that neutrinos 
should have non-zero masses to resolve this paradox. Soon after that work several 
authors ^, |^ made claims in contradiction to this conclusion. Smirnov and 
Vissani ||^ argued that the mechanism by which the neutrino sea in neutron stars 
leading to blocking of long-range forces due to the Pauli principle should resolve 
the paradox. Abada et al. argued that such self-energy is small, exactly solving 
a (1+1) dimensional potential model and a (3+1) dimensional fiat boarder model. 
Kiers and Tytgat made numerical analysis of the self-energy and the neutrino 
number of the ground state of a spherical neutron star, and found they tend to the 
semi-classical values 



W., ^ (1.1) 

^ '-^ (1.2) 

in the large volume limit and hence that no paradox exists. 

In this paper we calculate the self-energy and the neutrino number of a spherical 
neutron star without recourse to numerical analysis and derive the same conclusion 
as Kiers et al. based on a large volume approximation. In Section 2, we introduce the 
Schwinger formula for the Weyl spinor and explain our approximation to calculate 
the self-energy and the neutrino number in the neutron star. In Sectin 3, we give 
the results of our calculation. 



2. Formulation of Energy and Neutrino Number 

Let us first introduce the Schwinger formula for the ground state energy, W , in 
the case of a two-component (Weyl) spinor system. This is given by the difference 
between the ground state neutrino energy of the "vacuum" containing a neutron 
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star, |0), and that of the true vacuum |0). In the former, the neutrino propagates 
in the vector potential inside the neutron star due to exchange, 

l^=(0|//|0)-(0|//o|0), (2.1) 

where Hq is the free Hamiltonian, and H is the total Hamiltonian for the neu- 
trinos with the neutron star. The latter corresponds to the Lagrangian for the 
(2-component) neutrino field x and its hermitian conjugate x, 

L = xio^d,x^^Xo'X3^. (2.2) 

with cr'' = ctq = 1, and cr* = — cr,, (z = 1, 2, 3) the Pauli matrices. Here is the weak 
current of the neutrons proportional to the neutron density, 

~ {n^n)g^,Q. (2.3) 

In the following we consider the model Lagrangian of the neutron stars [^] 

£ = x(^a^a^ + M^))x, (2.4) 

where 

ix{x)= iie[R-\x\), (2.5) 

with R the neutron star radius and /x = ^^{ri)n) a constant, typically on the order 
of several eV. 

d 

Since the Hamiltonian densities in Eq. ( p.l| ) are given by — i lim — — trx(a;)x(2;'), 
the Schwinger formula for the energy is given by 

W = -i [ d^x^tT(S(x,x') - S°(x,x')).^^^^ (2.6) 

J OXq 

= -^JdETT{\n{E-H)-\n{E~Ho)), (2.7) 

where the symbol Tr is the trace over both spinor and configuration space indices, 
S{x, x') is the Feynman propagator defined as 

{iad + ii{x))S{x, x) = i6{x, x'), (2.8) 

{0\Tx{x)x{x')\0) = S{x,x'), (2.9) 

and S'°(x,x') is that of the free neutrino. In Eq. ( [^.T] ) and thereafter we use the 
same notation H for the quantum mechanical Hamiltonian, H = —icrV — fi{x), as 
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for the field theoretical one in Eq. ( |2.1| ), since this should have no confusion. To 
derive Eq. (|2.7|) from Eq. (|2.6|) we have used the formula 



S X ^ 



If 1 
— lim / dE(x\ — -, 



\x')e' 



-iE{xO-x'0) 



(2.10) 



and a similar one for iS'°(x, x'). The prescription with ie is introduced here to satisfy 
the boundary condition of the time ordering propagator. The neutrino number, q, 
is similarly given as 



q = — I (fxtT{S{x,x') — S'^{x,x'))x' 

±fdE'Tr(^ 

2-KiJ \E-H E-Hq 



(2.11) 
(2.12) 



In order to avoid the apparent divergence of the expressions for W and q, we 
differentiate the energy with respect to fi twice and the neutrino number once to 
obtain 



dq 1 



E-H 



y,x 



0{R-\y\), 



dfb 



2m 



dE / d^x / d^y 



E — H J x,y \E — H / y x 



9{R- \x\). 



(2.13) 
(2.14) 



In Eq. ( p.l3|) we see that the neutrino propagates from a point x inside the 
neutron star to another point y also inside the neutron star. Noting fiR ~ 10^^, we 
see that the neutrino propagator rapidly oscillates in the neutron star, and thus only 
a short distance with |a; — y| -C -R or finite n\x — y\ contributes to the integration of 
Eq. ( p.l3|) and (|2.14]) . Thus it is safe to replace the full propagator of the neutrino 



E - Ho + fie{R-r) 



x,y 



with the approximate one 



(2.15) 



(2.16) 



\E - Hq + hJ 

to evaluate Eqs. ( |2.13D or ( |2.14D , neglecting the small contribution of such neutrinos 
that virtually propagate out of the neutron star (see Fig. 1). 
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Figure 1: We neglect the neutrino propagating through path (b). We only take the 
path (a) into account. 

This approximation^ leads to 



E - Hq + ii9{R-r) 



where r = x - y and p+ = e{E){E + //) {e{E) = 1 (-1) for ^ > < 0)). The 
Eq. ( p.lTp represents nothing but a full propagator in condensed matter of infinite 
volume. 

The Schwinger formula (|2.7| ) has an ultraviolet divergence proportional to 
which corresponds to just the familiar vacuum polarization diagram for Z° and 
should be subtracted.^ Thus, the right-hand side of Eq. ( p.l3|) , a second derivative 
of W, has a constant divergence which should be subtracted. Equation ( |2.12| ) for the 
neutrino number also has a divergence. This divergence, linear in fx, corresponds to a 
one-loop correction to the neutrino density operator, familiar in the correction in the 
operator product expansion, and is caused by the operator mixing of the neutrino 
density and neutron density. This is also to be subtracted by renormalization. Thus 
the right-hand side of Eq. ( |2.14| ), the first derivative of q, should be renormalized to 
be zero at /i = 0. 



E-Ho + fi^ 



x,y 



(27r)3 E + ^-crp 



ip{x-y) 



^ err , . 

E + fi + — - 

r \r 



* Here the ie prescription is similar to that used in Eq. ( 2.10 ). 

^ We note that a divergence cannot be subtracted in case of an electron-condensate medium, 
since the exchange of charged weak boson is involved. In such a case we need a cut-off scale Ac, 
which is the weak boson mass or the inverse of the mean distance of the electrons. 
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3. Evaluation of Energy and Neutrino Number 



In this section, we evaluate the self-energy and neutrino number of a neutron 
star with Eqs. (|2.13|) and (|2.14|) using the approximation Eq. (|2.17|) . It is easy to 
formally obtain 



(fW 



d^x d?y r ^(cos 2/ir + /xr sin 2//r), 

loTT'^J \X\<R J \y\<R 



(3.1) 



with r = y — X. We renormalize it by subtracting the divergent term to obtain 



dfj.'^ 



d'^W d^W 



djj? 
1 



dfj? 
d^x 



(3.2) 



d^x d^y r ^(cos 2/ir + /ir sin 2/ir — 1). (3.3) 

Stt-^J \x\<r J \y\<R 



Noting the convergence of the integrand of Eq. (|3.3|) both at r = and r = oo, we see 
that the change of the integration region / ^xikr*^^^ I \y\<R'^^^ ~^ I \x\<Rd^^ I |r|<oo^^^ 
is allowed in the large volume limit. Then the integration in Eq. ( p.3|) is analytically 
performed to give 



d^W 



d^j? 



87r3 
1 



/ d^x I Aur^dr r ^(cos 2^r + jj,r sin 2yur — 1) 

J \X\<R J r<oo 



d^x{-4:7Ti/ 



Sir^J \x\<R 
2 

By integration we obtain the renormalized self-energy 

1 



W 



4 d3 



IStt 



(3.4) 
(3.5) 
(3.6) 

(3.7) 



where we have used the fact 
d^W 



Eq. (0)- 



that W is an even function of fi and VF|;x=o = 
0. The result Eq. ( |3.7| ) coincides with the semi-classical value of 



Here we make some remarks regarding some subtleties of the calculation. If we 
differentiate W with respect to fi more times, we obtain better convergence at short 
distance. In this sense it is safer to evaluate W in terms of higher derivatives. In 
fact we have no divergence in W^(")(/i) for n > 3. In the case of W^^\fi), apparently 
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we may have divergence worse than in the case of W^'^\ When we calculate W^^\fi) 
and use the approximation Eq. ( p.l7|) , however, we obtain a finite result, 



dW 



— dE d^x2(E + a) L^o (3.8) 

m J-oo J\x\<R 47rr 



dix 27ri J-oo J\x\<R 



' f d^x f^rcos^r -sin^r^^^^ ^^ ^^ 



2n'^ J\x\<R 

-^fi'R', (3.10) 



which is consistent with the renormalized one, Eq. (|3.7] ) . The ultraviolet divergence 
linear in // is accidentally cancelled in this calculation. Since Eq. (|3.9| ) has a superfi- 
cially divergent dimension, it requires a more careful calculation by taking account 

d'^W 

of the boarder effect. In order to study the boarder effect, let us take - — — , 

dfjdR 

d^W 1 r / 1 



/^^Tr(-^— ^ 5iR-\x\ 



+ I dE ti I d'x I d'yi——) 6{R-\y\) [——] 9{R-\x\). 



dfxdR 27rz 

2m J J " 'J \E-HJx,y"'~ \E-Hjy,x 

(3.11) 

We note this expression is exact. We again approximate it with the full propagator 
by Eq. ( [2.171 ). The result has an ultraviolet divergence linear in /i. 

After renormalizing as explained above, we obtain 

d-'W 2 , , cosAfiR + AfiR SijAfiR) 

iHdR "3^^ ^ + ^ 8^^:r ' ^^-^^^ 

where Si(x) = / dt . This gives the self-energy of the neutron star as 

Jo t 

W^-^4^{,^0m, (3.13) 



IStt V \R 

which again coincides with the classical value Eq. ( |1 . 1| ) in the large R limit. 

dQ 

Now let us evaluate the neutrino number. Equation (p.l4|) for — resembles 

Ct/i 

d^W 

Eq. ( |2.13| ) for ^ ^ , except for the sign and the range of the y integration. Since 
the quantity 

d'xf d^yiJ^—] (-^] (3.14) 



|a;|<R J\y\>R \E - H J x,y \E - H J y^x 



should be negligibly small for large fiR, they should coincide (except for the sign) 
to give 



dq 
dfi 



Sir 



with the aid of Eq. (|3.6| ). This gives the neutrino number, 

q = 

9n 



(3.15) 



(3.16) 



which is equal to the semi-classical value in Eq. (1.2 



4. Discussion 



In this paper we estimated the self energy and the neutrino number of neutron 
stars without recourse to numerical analysis, and confirmed that they coincide with 
the semi-classical values in the large volume limit. They do not lead to paradox 
nor to any lower bound for the neutrino masses. While our estimation was made 
non-perturbatively, it may be instructive to evaluate the self-energy perturbatively 
by expanding the integrand of Eq. (p.3|) with respect to /i and integrating each term 
over X and y without using the approximation Eq. ( ^.4[ ). The integration is easy 
to perform and leads to an alternating series in powers of fiR. This series is similar 
(though not exactly equal) to that considered in Ref. |^ and an example where only 
the sum (but not each term individually) is meaningful in large R limit. 



Note Added in Proof 



In Ref. Abada et al. estimated the self energy and the neutrino number 
by independent methods and got results different from ours. Their self energy is 
linear in the external potential (/i in our notation), while it should be even to keep 
CP invariance. We think that the difference is due to their incomplete Hamiltonian 
violating CP invariance without taking the average of it and its CP conjugation. Our 
method keeps the invariance, though implicitly, by taking the neutrino propagator 
with adequate ie procedure. Kiers et al. ^ also analyzed them independently with 
the result in good agreement with ours. 
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